Phase Diagram of mixed bond Ising systems by use of Monte Carlo and the 

effective-field theory 



O 

(N 



I 

C/3 



c3 



i 

c 

o 
o 



> 
o 
o 



(N 

o 
o 



J. B. Santos-Filho 1 , N. 0. Moreno 2 , Douglas F. de Albuquerque 3 

1 Departamento de Fisica, Universidade Federal de Minas Gerais, 31270-901, Belo Horizonte, MG, Brasil 
2 Departarnento de Fisica, Universidade Federal de Sergipe, 49100-000, Sao Cristovdo, SE, Brasil and 



Departamento de Matemdtica, 



Universidade Federal de Sergipe, 
(Dated: February 2, 2010) 



49100-000, Sao Cristovdo, SE, Brasil 



The phase transition of a random mixed-bond Ising ferromagnet on a cubic lattice model is 
studied both numerically and analytically. In this work, we use the Cluster algorithms of Wolff and 
Glauber to simulate the dynamics of the system. We obtained the thermodynamic quantities such 
as magnetization, susceptibility, and specific heat. Our results were compared with those obtained 
using a new technique in effective field theory that employs similar probability distribution within 
the framework of two-site clusters. 
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I. INTRODUCTION 

The study of the effects of disorder in magnetic systems 
has been an object of intense investigations during the 
last five decades. 

The Monte Carlo technique is an useful tool which in 
many cases, gives better results regarding other meth- 
ods from analytical approximations. The influence of 
quenched, random disorder on phase transitions is of 
great importance in a large variety of fields [l| . For pure 
systems exhibiting a continuous phase transition, Har- 
ris @ derived the criterion that random disorder is a 
relevant perturbation when the exponent of the specific 
heat of the pure system is positive, a > 0. In this case 
one expects that the system falls into a new universality 
class with critical exponents governed by a disordered 
fixed point. For a < disorder is irrelevant, and in the 
marginal case a = 0no prediction can be made. 

Since for the three-dimensional (3D) Ising model it 
is well known that a > 0, quenched, random disorder 
should be relevant for this model. In three dimensions 
(3D) most of the computer simulation studies have con- 
centrated mainly on the site-diluted Ising model 0, [HJ . 

In this work, we study the Ising model with mixed- 
bond by using of Monte Carlo simulation, applying the 
algorithm cluster of Wolff Q . 



II. MODEL AND SIMULATION SETUP 

We study the spin 1/2 ferromagnetic Ising mixed-bond 
model defined by the following Hamiltonian 



(ij) 



fa = ±1), 



(1) 



where the sum extends over all pairs of neighboring sites 
on a cubic lattice of linear size L with periodic bound- 
ary conditions, ft — l/fegT and the exchange couplings 



Kij are allowed to take two different values 



K 



J/ksT and 0. The interactions are assumed to be inde- 
pendent random variables with distribution 



P(K i:j ) = P S(K tJ -K) + {1- p)8{K lJ - XK), (2) 

where p is the concentration of magnetic bonds in the 
system bonds such that p = 1 corresponds to the pure 
case and A is the competition parameter with |A| < 1. 

The simulations were performed on a set of following 
lattice sizes L = 10, 16, 20, 26, 30, 36, 40 with peri- 
odic boundary conditions. The aim of the first set of 
simulations is to estimate the critical temperature of the 
model at different L. Due to the finite-size scaling theory 
@ , the finite system of linear size L will demonstrate an 
evidence of a critical behavior at a certain temperature 
Tc{L) which differs from the critical temperature of the 
infinite system Tc(oo) Q. 



T C {L) =T c {oo) + OiL' 1 ''" + 



(3) 



where the correction-to-scaling terms have been omitted. 

The static thermodynamic quantities of interest in- 
clude the average magnetization M and the magnetic 
susceptibility \ 



X 



[<M 2 ) - (M) 5 



(4) 
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The phase diagram is obtained numerically from the 
maxima of a diverging quantity. Here we choose the mag- 
netic susceptibility, since the stability of the disordered 
fixed point implies that the specific heat exponent is neg- 
ative in the random system pi, 9]. Thus, the error in this 



2 



quantity is larger than for the susceptibility. To get an ac- 
curate determination of the maxima of the susceptibility, 
we used the histogram reweighting technique with 2500 
Monte Carlo sweeps (MCS) and between 2500 and 5000 
samples of disorder. The number of Monte Carlo sweeps 
is justified by the increasing behavior of the energy auto- 
correlation time, te, and we chose for each size at least 
250 independent measurements of the physical quanti- 
ties (Nmcs > 250 te)- The choice of Nmcs is justified 
by the increasing behavior of the energy autocorrelation 
time te as a function of p and L. At the critical point of 
a second-order phase transition one expects a finite-size 
scaling (FSS) behavior te oc L z , where z is the dynamical 
critical exponent 



X= -0.4 




III. RESULTS 

The figure Q] display curves of magnetization versus 
temperature through computational simulation (Wolff) 
for p = 1 and L = 15, 20 and 30. The critical temperature 
obtained when p = 1 was Tq = 4.510 it is close of the 
expected value Tc = 4.51. [T3| We observe that the curves 
keep the same behavior, in spite of in the proximities 
of the critical point they move away each other. The 
critical point was estimated of the inflection of the curve. 
It can be notice that increasing the lattice size, i.e. the 
value of L, we get more precision to estimate the critical 
temperature. 

The magnetic susceptibility as a function of temper- 
ature for different p values and for A = —0.4 and L = 
40 is shown in Fig. [2J The peaks are sharper for values 
lower of p. To p < 0.65 and A = —0.4 the susceptibility 
not displays the peak associate with the magnetic tran- 
sition due to competitions of the exchange interactions. 
When A > the system presents always long range order. 
We used the histogram reweighting technique with 2500 
Monte Carlo sweeps (MCS) and between 2500 and 5000 
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FIG. 1: Magnetization vs. ksT / J p = 1 for the Ising model 
with L = 15, 20 and 30 in a cubic lattice. 



FIG. 2: x vs - ksT/J for the Ising model with A 
L = 40 for several concentrations p. 



-0.4 and 



samples of disorder to get an accurate determination of 
the maxima of the susceptibility. 

The phase diagram obtained from the location of the 
maxima of the susceptibility for the largest lattice size 
(L = 40) as a function of the concentration of magnetic 
bonds is shown in Fig. [3] for A=0, 0.1, 0.2, 0.5, -0.2 and 
-0.4. Solid lines are the predictions of the effective field 
approximation. A very good agreement with the simu- 
lated transition line is obtained. 
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FIG. 3: Phase diagram of the 3D mixed-bond Ising model 
compared with the effective field approximation. 

The reduced fourth-order Binder cummulants [j~lT| sup- 
ply an alternative method to estimated critical points 
that can be determined from the crossing point of the 
cumulants for different lattice sizes. It is calculated by 



U L = 1 



(m 2 y 



where [...] denotes the average over disorder and (...) 
refers at the thermal average. As an example, we show 
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FIG. 4: Fourth-order Binder cumulant Ul vs. fesT/J, for 
different lattice sizes as indicated in figure. 



in Fig |4] the T- dependence of the reduced fourth-order 
Binder cummulant for p = 1 and for various lattice sizes. 
Critical temperature obtained from this figure is in agree- 
ment with those obtained of the maxima of the magnetic 
susceptibility. 

The peak locations of the maxima susceptibility for 
each L are plotted versus L" 1 ^, where the value of v 
is determined of the linear fit of log-log plot ^Jpjr vs L 
(not shown here). The critical temperature can be es- 
timated from an infinite-size extrapolation in according 
with Eq.(3). We illustrate this procedure in Fig [5] for 
A = 0.5 and p = 0.4, the fit yield to T c = 3.1065(3). 



TABLE I: Critical exponents for three dilutions and A = 0.5. 



p 


1 


0.8 


0.6 


0.4 


l/v 


1.591(7) 


1.68(7) 


1.63(4) 


1.56(3) 


P/v 


0.4920(5) 


0.49(1) 


0.50(13) 


0.496(5) 




2.008(6) 


2.06(3) 


1.950(6) 


2.001(6) 



3.106 



3.104 



O 



3.102 



p = 0.4 

X = 0.5 



The average magnetization m and the magnetic sus- 
ceptibility x scale with the lattice size as: 



rn 



,L 



-p/v 



(6) 



here a m and a x are non-universal amplitudes. From these 
power-laws we extracted the exponents (3 and 7 plotting 
in logarithm scale the lattice size dependence of the sus- 
ceptibility and average magnetization. The critical expo- 
nents obtained to p = 0.4, 0.6, 0.8 and A = 0.5 are shown 
in Table I. The critical exponents oscillate without hav- 
ing an apparent correlation with the dilution. They are 
pretty close of those for disordered Ising model [l2[ 



IV. CONCLUSION 



We carried out Monte Carlo simulations for study the 
influence of bond dilution on the critical properties of 
the Ising Model applied for cubic lattice. We obtained 
thermodynamic parameters for |A| < 1. 

Satisfactory results are obtained using the algorithm of 
Wolff and showed that this technique is appropriated to 
treat the Mixed-bond problem. The Monte Carlo tech- 
nique results gives similar results to the obtained ones by 
the effective field theory. 

The critical behavior of the mixed-bond model is gov- 
erned by the same universality class as the site-diluted 
model and pure Ising model. 



3.100 



0.00 



0.01 



0.02 
-1/v 



0.03 



ACKNOWLEDGMENT 



FIG. 5: Size-dependent critical temperature Tc(L) vs. L 1/7 " 
for A = 0.5 and p = 0.4. 



N.O. Moreno is partially supported by FAPITEC-SE 
and CNPq Brazilian agencies. J. B. Santos Filho ac- 
knowledge financial support by CAPES. 



4 



[1] R. Folk, Y. Holovatch, T. Yavorskii: Physics Uspiekhi 

173 (2003) 175. 
[2] A. B. Harris: J. Phys. C 7(1974) 1671. 
[3] S. Wiseman and E. Domany, Phys. Rev. Lett. 81, (1998) 

22. 

[4] H.G. Ballesteros, L.A. Fernandez, V. Martin-Mayor, A. 

Mufioz Sudupe, G. Parisi, J.J. Ruiz-Lorenzo, Phys. Rev. 

B 61, (2000) 3215. 
[5] U. Wolff, Phys. Rev. Lett. 62 (1989) 361. 
[6] M. N. Barber Finite-size scaling. In: Phase Transitions 

and Critical Phenomena, eds. C.Domb and J.Lebowitz, 

Academic, New York, 1983, vol.8 
[7] A. M. Ferrenberg, D.P Landau, Phys. Rev. B 44 (1991) 

5081. 



[8] J. T. Chayes, L. Hayes, D. S. Fisher, T. Spencer, Phys. 

Rev. Lett. 57 (1986) 2999. 
[9] L. E. Zamora, G. A. Perez Alcazar, A. W. Bohorquez, 

A. A. Rivera and J. A. Plascak, Phys. Rev. B 51 (1995) 

9329. 

[10] C Domb and M. S. Green, 1976. Phase Transitions and 

Critical Phenomena, ed. por C. Domb e M. S. Green, 

(Academic, London) 
[11] K. Binder and D. W. Heermann, Monte Carlo Simulation 

in Statistical Physics. An Introduction, 3rd ed. (Springer, 

Berlin, 1997). 

[12] H. G. Ballesteros, et al, Phys. Rev. B 58 (1998) 2740. 



